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A CHANNEL ESTIMATOR 

CROSS-REFERENCE 

[001] This application is a Continuation Application of U.S. 
Patent Application Serial Number 09/438,475, filed on November 12, 1999. 

FIELD OF THE INVENTION 

[002] The present invention relates to channel estimators in wireless communication 
systems generally and to such channel estimators that operate on data and pilot signals in 
particular. 

BACKGROUND OF THE INVENTION 

[003] In practical digital communication systems, the frequency response of the 
underlying channel from the transmitter to the receiver is seldom known at the receiver 
side. For example, in digital communication over the dial-up telephone network, the 
communication channel will be different every time a number is dialed, because the 
channel route will be different. In this example, the characteristics of the channel are 
unknown a priori. There are other types of channels, e.g. wireless channel such as radio 
channels and underwater acoustic channels, whose frequency response characteristics 
are time varying. Thus coherent communications for such channels require the 
utilization of adaptive algorithms, known as "channel estimators", for tracking / 
estimating the varying characteristics of the channel. 

[004] Traditionally, channel estimators are divided into two categories: data aided 
and non-data aided (blind) estimators. Data aided channel estimators operate on a pre- 



1 



Attorney Docket Number: P-2654-US1 Assignee: Intel Corporation 

specified set of transmitted symbols that are known to the receiver. These symbols do 

not convey any information and are often called cc pilot symbols" or 'Gaining sequences". 

Data aided channel estimators are typically simple to implement and relatively robust. 

Their major disadvantage is that they lead to an overall reduction in system throughput, 

since some of the transmitted symbols (the pilot symbols) do not carry any information. 

[005] Non data aided channel estimators, on the other hand, do not reduce the 

system throughput. However, they are typically quite complicated to implement as they 

are often based on higher order moments/cumulants of the received signal, and they 

most often suffer from high statistical variability, i.e. they suffer from large estimation 

errors. 

[006] The article "Maximum A Posteriori Multipath Fading Channel Estimation for 
CDMA Systems" by Mohamed Siala and Daniel Duponteil, Proceedings of Vehicular 
Technology Conference. Houston, Texas, May, 1999, describes a channel estimation 
algorithm which combines both approaches. This algorithm uses both pilot and data 
symbols to construct a channel estimator. However, this algorithm requires that the joint 
statistical probability distribution of the channel multipaths be known to the receiver. In 
practice, a complete statistical description of the channel characteristics is seldom 
known to the receiver. Moreover, these characteristics may be time varying. 
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SUMMARY OF THE INVENTION 

[007] An object of the present invention is to provide an improved channel 
estimator without using any a priori statistical information about the channel. Instead, 
the present invention uses a priori probabilities of the received symbols, be they pilot, 
data, power control, etc. symbols. 

[008] There is therefore provided, in accordance with a preferred embodiment of 
the present invention, a channel estimator based on the values of received data and on a 
priori probabilities only of received symbols. The channel estimator includes a symbol 
probability generator, a noise variance estimator and a channel tap estimator. The 
symbol probability generator generates a priori probabilities only of transmitted symbols 
found in the received signal(s). The noise variance estimator estimates at least one noise 
variance corrupting the received signal(s). The channel tap estimator generates channel 
estimates from the received signal(s), the a priori probabilities and the noise variance(s). 
[009] Additionally, in accordance with a preferred embodiment of the present 
invention, the channel tap estimator solves the following equation: 

hML=^'Sy(t)-z(t;hJ 

[0010] Moreover, in accordance with a preferred embodiment of the present 
invention, the channel tap estimator includes a z-unit, a combiner and a channel tap unit. 

The z-unit generates a z-value for zteh^) from the a priori probabilities, the noise 
variance(s) and the channel estimates. The combiner combines the z-value with the 
received signal(s). The channel tap unit determines channel tap values from the output 
of the combiner. 
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[001 1] Further, in accordance with a preferred embodiment of the present invention, 

the z-unit calculates the following equation for a quadrature phase shift keying (QPSK) 

channel: 



z(t;h)» 



alx/i *\ 4Real{<7(tXh->(H)} - ^Real{<y(t),h->(l + j)} 

Pi(t)(l + J)e° +P 2 (t)(l-J)e a 



-. 4-Real{<F(t),h->(-l-j)> -lReal{<y(t),h->(-l + j)} 

+ P 3 OX" 1 + J) e ° + P4 (t)(~l - J)e° 



4-R«al{<y(0,h->(l-j)} , . 4"R«a>{<y(0,h , >(l+j)} 

Pi(t)e° +P 2 (t)e° 



, , ^Real{<y(t),h->(-l-j)} 4-Real{<y(0,h-X-l+j)}^ ' 

+ p 3 (t)e° +P 4 (t)e° 



T 



where y (t) is a vector of the received signal(s), h is a vector of the channel 
tap estimates and p(t) is a vector of the symbol probabilities. 

[0012] Still further, in accordance with a preferred embodiment of the present 

invention, the channel tap unit is a summer over a window of received symbols. 

[0013] Additionally, in accordance with a preferred embodiment of the present 

invention, when s(t) is a pilot symbol, then p(t) = where the ^'s are the 

normal basis vectors having a one (1) in their i-th component and zero otherwise. When 

s(t) is a data symbol, then 

P, (t) = P 2 (t) = p 3 (t) = p 4 (t) = 0.25 

When s(t) is a transmit power control (TPC) symbols, then 

p 1 (t) = p 4 (t) = 0.5, p 2 (t) = p 3 (t) = 0 
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[0014] Moreover, in accordance with a preferred embodiment of the present 

invention, the channel tap unit includes an anchor unit, an averager and an interpolator. 
The anchor unit determines a pilot anchor h p using N p pilot symbols of one time slot 

n and a data anchor h s using N s data symbols of the time slot n . The averager 
averages the pilot and data anchors to produce a slot anchor ^anchor 

(n) and the 

interpolator interpolates between adjacent anchors n anc j lor (n-l) anchor 

(n) 

to obtain channel estimates for the II -th slot 

[0015] Further, in accordance with a preferred embodiment of the present invention, 
N s = 2N p and N p of the data symbols can be taken from before the pilot symbols 

of the time slot and N p of the data symbols can be taken from after the pilot symbols. 

[0016] Still further, in accordance with a preferred embodiment of the present 
invention, the unit for linearly interpolates includes unit for separately interpolates 
amplitudes and phases of the adjacent anchors. 

[0017] Additionally, in accordance with a preferred embodiment of the present 
invention, the channel estimator also includes a unit that averages slot anchors as 
follows: 

M /v 

Sh mch „ r (n-l-k)p k 
H mC h 0 r ( n_ l) - 



anchor V / m 

k=-M 



where H anchor is the average value of slot anchors h anchor in a slot ranging 

between the values -M to M for k, and f3 k is a user defined weight factor. 
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[0018] Moreover, in accordance with a preferred embodiment of the present 

invention, the z-unit includes a lookup table unit. 

[0019] Further, in accordance with a preferred embodiment of the present invention, 
the at least one received signal is the output of at least one despreader. Alternatively, it 
could be the downconverted and demodulated output of more than one antenna. 
[0020] Still further, in accordance with a preferred embodiment of the present 

2 

invention, the noise variance estimator generates a noise variance (T as follows: 

6 k 2 (n) = (1 - a) • 6 k 2 (n - 1) + a a k 2 (n) 

where the time index II is in units of slot, OC is a user selectable exponential 
forgetting factor, and 

2 1 N " 

N p t=i 

where N p is the number of pilot symbols per slot, G(t) is an automatic gain 

control (AGC) level and y k (t) is the k -th element of the vectory(t) . 

[0021] Further, in accordance with a preferred embodiment of the present invention, 
the channel tap unit is an infinite impulse response (OR) filter given by: 

h(t) = -h(t-1)a 1 -,...,-h(t-p)a p +y(t)z(t;h(t-1))*b 0 

+ y(t-1)z(t-1;h(t-2))*b 1 + y(t-q).z(t-q;h(t-q-1))*b q 

where SL i and bj are user-defined parameters as are the filter orders p and 

q 

[0022] Alternatively, there is also provided a channel estimator having two 
estimators, a first estimator that operates on a continuous pilot channel and a second 
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G(t) 



s(t)' 
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estimator that operates on a traffic channel with interleaved pilot symbols. The 

estimator also includes a combiner that combines the output of the first and second 

estimators. The first estimator is similar to those of the prior art while the second 

estimator is similar to that described hereinabove. 

[0023] Finally, the present invention includes the methods performed by the channel 
estimators disclosed herein. 
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BRIEF DESCRIPTION OF THE DRAWINGS AND APPENDIX 

[0024] The present invention will be understood and appreciated more fully from the 
following detailed description taken in conjunction with the appended drawings and 
Appendix in which: 

[0025] Figs. 1 A, IB and 1C are scatter plot illustrations of the values of a function, in 
the presence of varying levels of noise, where the signal to noise ratio (SNR) is high, 
medium and low, respectively; 

[0026] Fig. 2 is a schematic illustration of a batch mode channel estimator, 
constructed and operative in accordance with a preferred embodiment of the present 
invention; 

[0027] Figs. 3 A and 3B are schematic illustrations of implementations of the channel 
estimator of Fig. 2, for single and multiple CDMA fingers, respectively; 
[0028] Fig. 4 is a schematic illustration of an adaptive channel estimator, constructed 
and operative in accordance with a preferred embodiment of the present invention; 
[0029] Fig. 5 is a schematic illustration of two anchors, useful in understanding an 
anchor embodiment of the channel estimator of Fig. 4; 

[0030] Fig. 6 is a schematic illustration of an implementation of the channel 
estimator of Fig. 4, for a single finger; and 

[0031] Appendix A provides a series of approximations to a main component of the 
channel estimator to simplify its implementation. 
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DETAILED DESCRIPTION OF THE PRESENT INVENTION 

[0032] The present invention is a channel estimator for linearly modulated, digital 
communication schemes. It takes as input the downconverted and demodulated received 



despreaders in a code division, multiple access (CDMA) rake receiver (as known in the 

art) or the demodulator output in a PAM/QAM (pulse amplitude modulation / 

quadrature amplitude modulation) receiver (as also known in the art). 

[0033] The present text will emphasize CDMA cellular systems; however, the 

present invention can be applied to any digital, linear modulation scheme. 

[0034] It will be appreciated that, in a receiver having more than one antenna, 

different elements of the vector y(t) will correspond to the downconverted and 

demodulated signals from the different antennas. In a multiple antenna, CDMA 

receiver, there are multiple despreaders per antenna. The elements of the vector y (t) 

correspond to the output of all of the despreaders. 

[0035] Various embodiments of channel estimator 22 are shown in Figs. 2, 3A, 3B, 
4, 5 and 6. To understand the structure of the channel estimator of the present invention, 
a derivation is first provided. 
Derivation 

[0036] The vector output y (t) at time instant t is modeled by a simplified static 



channel model as a function of the transmitted signal S(t) , a vector of static channel 



signal, denoted by a vector 




For example, y(t) 



can be the output of 




taps h , and a vector n(t) of noise samples. The model is: 



y(t) = hs(t) + n(t) 
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Equation 1 

[0037] It will be appreciated that: 

dim{y(t)} = dim{h} = dim{n(t)} > dim{s(t)} = 1 

where dim{v} staads for the dimension of the vector V . Thus, the number of 
unknown channel taps is denoted by dim{h} which is at least 1. For example, in a 
CDMA rake receiver, dim{h} equals the number of active fingers. The noise term 
Il(t) for a CDMA system represents the additive white Gaussian noise (AWGN), the 
cross-talks from other users, cross-talks from un-modeled channel taps (i.e. channel taps 

that do not correspond to any of the active fingers), etc. and the channel estimate ll is 
the unknown finite dimensional parameter vector that needs to be estimated based on 

observation of y (t) for t = 1,2, . . . , T . 

[0038] In accordance with a preferred embodiment of the present invention, the 
following assumptions are made in order to generate the channel estimator: 

1 . S(t) is a sequence of independent and identically distributed (Li.d.) random 
variables of data symbols with time varying, a priori symbol probabilities 
p(t). For simplicity, S(t) will be considered herein as a quadrature phase 
shift keying (QPSK) signal. 

2. Il(t) is an i.i.d. sequence of complex valued, zero mean Gaussian random 

2 

vectors with independent components whose variance is CJ . Conventional 
circular Gaussian probability distribution functions (pdf s) are assumed, i.e. the 
real and imaginary parts of each of the noise components are independent with 
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variance CJ 2 / 2 . The case of unequal noise variances (that is, when each 
component of n(t) possesses different variance) is also considered 
hereinbelow. 

3. S(t) and Il(t) are mutually independent 
[0039] The present invention defines a probability vector p(t) that defines the a 

priori probabilities for the various symbols of s(t) . In CDMA systems, signals are sent 

between a mobile and a base station in pre-defined time slots. In each time slot, many 
symbols are sent, only some of which are raw data symbols, e.g. the compressed voice 
data. The remaining symbols are used for controlling the communication, estimating the 
unknown channel taps, etc. Among these control symbols there are the pilot symbols, 
which, as described hereinabove, are a set of symbols that are a-priori known to the 
receiver. The pilot symbols are used to synchronize the base and mobile stations and for 
estimating the channel taps. Other control symbols include transmit power control 
(TPC) symbols. With the current standardization regarding 3 rd generation cellular 
systems (3GPP RAN wideband CDMA), these TPC symbols are known a-priori to have 
only one out of 2 possible values. The TPC symbols are defined in the 3GPP document, 

RAN 25.214 VI. 1.1(1999-09). The probability vector p(t) allows the present 

invention to treat both pilot and data symbols. 

[0040] For a QPSK symbol constellation, p(t) is defined as follows: 
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Pr{s(t) = l+j} 




Pi(t) 


Pr{s(t) = l-j} 


= 


P2O) 


Pr{s(t) = -l + j} 




p 3 (t) 


_Pr{s(t) = -l-j}_ 




P4(t). 



Equation 2 

[0041] If s(t) is a pilot symbol, then there exist some 1, l<i<4 such that 

P(t) = e s 

Equation 3 

where the G^s are the normal basis vectors having a one (1) in their i-th 

component and zero otherwise. The specific vector used per time instant depends on 
the transmitted pilot symbol. For example, the following sequence of pilot symbols 

s(l) = 1 + j, s(2) = 1 - j, s(3) = 1 + j, s(4) = -1 - j 

Equation 4 

corresponds (from Equation 2) to the following sequence of e. 's 

p(l) = e, , p(2) = e 2 , p(3) = e, , p(4) = e 4 

Equation 5 

If, on the other hand, s(t) is an equally likely data symbol, then 

Pi(t) = P 2 (t) = p 3 (t) = p 4 (t) = 0.25 

Equation 6 

[0042] Furthermore, if there is some a-priori information on the transmitted data, it 
can be incorporated herein. For example, the transmit power control (TPC) symbols are 
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known a-priori to have only one out of 2 possible values either l+j or -1-j. For TPC 



symbols, the corresponding p(t) would be 



P(t) = 



0.5 

0 

0 

0.5 



Equation 7 

[0043] Now, with the above assumptions, the conditional pdf of y(t) given s(t) 



and h is Gaussian and it is given by, 



P(y(t)/s(t);h)=c.e 



|y(t)-h- S (t)f 



-VTly-cW-MCt)! 1 
= c • e kd 



Equation 8 



13 



Attorney Docket Number: P-2654-US1 



Assignee: Intel Corporation 



[0044] Therefore, 

P(y(t);h)=c 



= c 


Pl(t) 


1 Dlniiyi (1 

~ Z |y k (t)-h k -(i+j)r 

>e a M 


+ 


p 2 (t) 


1 DI Eiy»i ii 

~r £ y k (t)-h t (i-j) 

, e a w 


+ 


p 3 (t) 


• e c M 


+ 


p 4 (t) 


~^T , |y k (t)-h k -(-i-j)r 
. e CT 



Equation 9 

[0045] Due to the i.i.d nature of both s(t) and n(t) , the joint pdf of y(t) for 

t = 1,2,..., T is simply the product of the one dimensional pdf s, or in the log 
domain: 



T 

Loglp^hl^Log 



-V Z |y*(0-V(W)| 2 



-r Z bkW^h-O-yjf 

— x Z k(o-A* (-i+y)| l 
— t Z Mo-vc-W)! 1 



+ TLog[c] 



Equation 10 

where, for notational brevity, Y denotes the vector whose t-th sub-block is 
y(t) for t = 1,2,..., T (i.e. Y is a (Dim{y} • T) X 1 column vector). 
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[0046] Now, maximum likelihood (ML) estimation of h amounts to differentiating 

Equation 10 with respect to h and equating to zero. To simplify the notation, we define 

(J)(t) as follows. 



f 



+(t) = 



\ 



^=7 + J- 



^dReal{h} dlmag{h} 



Log[p(y(t);h)] 



[0047] 



[0048] 



Equation 1 1 

The ML channel estimator is obtained by solving the implicit equation: 

i?(t) = o 

t=l 

Equation 12 

Using the equality below, 

,^ -4 £ y k (t)-h k -0+j) 2 M -4 £ |y k (t)-h k d-j)| 2 

Pl (t)-e * ~ + p 2 (t)-e a ~ 

,^ -4- £ MO-M-i+j)' ~ 2- |y k (t)-h k (-i-j)r 

+ p 3 (t)e " + p 4 (t)-e 



= e ff 



yiy ( t)r + 2|h|f) 



Pi(t)-e« 



r Real{<y(t),h*>-(l-j)} 



4ReaI{<y(t),h*>.(l+j)} 

+p 2 (t).e° 

4-R ea H<y(t),E'>(-i-j)} 
+ P 3 (t)-e CT 

4R ea H<y(t),E*>(-i+j)} 
+ P 4 (t)-e CT 



Equation 13 

we obtain, after quite a few algebraic manipulations: 
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<t>(t) = -4r • h + ~ y(t) • z(t; h) - 
a a 

Equation 14 

where (•)* stands for the conjugate of the bracketed expression and the scalar 
process z(t;h) is given by: 



z(t;h)^ 



✓ w, «\ ^Real{<y(t),h>(l-j)} 

Pl (t)(l+j)e' +P 2 (t)(l-J)e 



. x 4Real{<y(t),h , > 



+ P 3 (t)(-l + j)e 



>n ^Real{<y(t),h'>(-l-j)} A Re al{<y(t) J 

+ p 4 (t)(-l-j)e ff 



X 



Pi(t)e 



4-Real{<y(t),h'>(l-j)} 4- R ^ al {<y(t),h*>(l+j)} 

+ p 2 (t)e CT 



-^-Real{<y(t),h*>(-l-j)} 4 ReaI {<y( t ).h*>(-l+j)} 

+ P 3 (t)e CT +P 4 (t)e° 



-1 



Equation 15 

< X, y > denotes the scalar product between the vectors X and y . 
Combining Equation 12 and Equation 14, the ML estimator becomes the solution 



of the implicit equation: 



hML=^-iy(t)z(t;h ML )* 



2T t~ 

Equation 16 

[0049] It will be appreciated that, in the general setting, no closed form solution 
exists to Equation 16. However, there are some special cases of the ML estimator of 
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Equation 16. First, consider the case of pilot symbols, by substituting Equation 3 into 

Equation 15: This results in: 

z(t;h) = s(t) 

Equation 17 

and Equation 16 reduces to: 

l x 



h ML = 



Ey(t)-s(t)* 



2T t~" 

Equation 18 

which is the conventional data-aided ML channel estimator. 

In the non-data-aided setting where Equation 6 holds, Equation 15 reduces to: 



z(t;h) = 



(1 + j) ■ sinhf 4-Real{< y(t),h* > (1 - J)} 



(2 



+ (l-j)-sinh — Real{<y(t),h*>(l + j)} 



x 



r 



cosh 



Real{<y(t),h*>(l-j)} 



J 



-i-i 



+ coshf4-Real{<y(t),h* >(l + j)}^ 

Equation 19 

[0050] In a high SNR setting, it can be shown that, in both nominator and 
denominator of z(t \ h) , only one of the exponents is dominant and hence, 

z(t;h)«s(t) 

Equation 20 
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where s(t) is the nearest neighbor to < y(t),h > from the source 

alphabet (i.e. §(t)is the hard decision of < y(t),h > ). In this case, the ML 
estimator reduces to: 

2T t=i 

Equation 21 

which is the well known, decision directed, channel estimator. 
[005 1] La a low SNR setting, the ML estimator of the present invention, in effect, uses a 
soft symbol metric that is based on the instantaneous SNR and follows a hyperbolic law. 
To illustrate this, reference is now made to Figs. 1A, IB and 1C which graph the 
function in Equation 19 for several SNR values (20dB, lOdB and OdB, respectively) and 
for the scalar case of one active finger. Each figure shows a scatter plot of the values of 

z(t \ h) in the complex plane. 

[0052] As can be seen, at the high SNR of Fig. 1A, z(tj h) always takes on one of 
the four hard decision values (1 - j), (1 + j), (-1 - j), (-1 + j) of the received signal. This 

matches with Equation 20 and 21. However, at lower SNR's, z(t;h)is "softened" 
according to the confidence of each specific hard decision. Thus, Fig. IB (medium 
SNR) has values for z(t \ h) mostly along the real and imaginary axes and Fig. 1C 

(low SNR) has values for z(t \ h) all over the complex plane. 

[0053] Until now, the noise variance has been identical for all elements of the noise 
vector H(t) . It is also possible to relax this assumption. Thus, in Equation 15 and 19, 
the following 
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4rReal{<y(t),F> } 

a 

Equation 22 

is replaced by 

2Real{<y(t),h*> } 

Equation 23 

where y (t) is obtained from y (t) by dividing each of its components by 
the corresponding noise variance, i.e. 

{ y(t) }l .M>k 
a i 

Equation 24 

Appendix A provides a series of approximations to Equation 19 to simplify its 
implementation. 

[0054] In some applications, such as the 3 GPP wideband CDMA cellular systems, 
there is a continuous pilot channel separate from a traffic channel that contains data and 
pilot symbols. For these applications, a channel estimator based on both the traffic and 
the pilot channels can be constructed. The channel estimator of the present invention can 
be used for the traffic channel while a prior art channel estimator can be used for the 
continuous pilot channel. The two channel estimates can be combined to produce the 
final channel estimate. This provides a statistically more stable channel estimator. 
Implementation 
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[0055] Reference is now made to Fig. 2, which schematically illustrates a channel 

estimator 30 of the present invention. Estimator 30 comprises a noise variance estimator 

32, an a priori symbol probability generator 34 and a channel tap determiner 36. 

2 

[0056] Noise variance estimator 32 determines the noise variance CT as described 
hereinbelow. There are a variety of methods, known in the art, for estimating the noise 
2 

variance CJ . One method, incorporated herein by reference, is described in the article 
"An Efficient Algorithm for Estimating the Signal-to-Interference Ratio in TDMA 
Cellular Systems", by Mustafa Turkboylari et al., IEEE Transactions on 
Communications. Vol. 46, No. 6, June 1998, pp. 728 - 731. 

[0057] Another method of estimating the noise variance, suitable for the 3GPP 
wideband CDMA time slot structure, uses pilot symbols, only and implements the 

k-^-2 
-th element of the noise variance vector CJ : 

6 k 2 ( n ) = (1 - a) • d k 2 (n - 1) + a • a k 2 (n) 

Equation 25 

[0058] Where Hie time index II is in units of slot, OC is a user selectable exponential 
forgetting factor, and 



2 1 N - 

5 k 2 (n)-—E 
N p t=i 



y*ct> s(t) - 



i &y k (t) 



2^. S (t)* 



N p <~ G(t) 



G(t) 

Equation 26 

where N p is the number of pilot symbols per slot, G(t) is the receiver's 
automatic gain control (AGC) level (as known in the art) and y k (t) is the k -th 
element of the vector y(t) . 
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[0059] Symbol probability generator 34 receives information from higher layers in 

the receiver defining the type of the current symbol. For example, the higher layer may 

indicate that the current symbol is a pilot, a power control or a data symbol. Symbol 

probability generator 34 then produces the probability vector p(t) associated with the 

symbol type, as per Equations 2-7. 

[0060] Channel tap determiner 36 determines the channel tap estimate vector by 
solving Equation 16. Since Equation 16 is an implicit equation, one has to resort to 
iterative algorithms for solving it. There are numerous iterative approaches, based on 
gradients and/or Hessians of Equation 16 that can be implemented. See, for example, 
the book Numerical Recipes in C: The Art of Scientific Computing, by Press et al., 2 nd 
Edition, Cambridge University Press, 1992. 

[0061] Channel estimator 30 operates in a "batch" mode, taking a sequence of T 
samples (where T is a user selectable parameter often related to the fading rate) and 

iterating until the channel estimate h is stable. 

[0062] As with any iterative algorithm, an initial point must be provided for the 
algorithm. There are a variety of initialization procedures. For example, one may 

arbitrarily set the initial h to unity and then apply the iterative algorithm to converge to 
the true value. Alternatively, one can apply a conventional channel estimation technique 

based only on pilot symbols and use its output as the initial value for h . 
[0063] Reference is now made to Fig. 3A, which schematically illustrates the 
implementation of a simple iterative solution to Equation 16 that does not require 
calculation of gradients or Hessians. The channel estimator 30 of Fig. 3 A is a single tap 
channel estimator. 
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[0064] In the embodiment of Fig. 3 A, channel tap determiner 36 comprises a scalar 

z-generator 37, a multiplier 38 and a summer 39 over a period of length T of the output 

of vector multiplier 38. 

[0065] Z-generator 37 receives the demodulated scalar signal y(t) , the probability 
vector p(t) from probability generator 34, the noise variance CT from noise variance 

estimator 32 and the channel estimate h' of the previous iteration. Z-generator 37 then 
generates the scalar value z(t;h') from Equation 15. For example, this operation can 

utilize a lookup table. In particular, when Equation 19 needs to be computed, a lookup 
table can be used which saves the need to calculate the hyperbolic sine and cosine 
functions. 

[0066] Multiplier 38 multiplies the output y(t) with the value z(t;h') and 

summer 39 sums the output of multiplier 38 over a period of length T. The result is the 

updated channel estimate h which is then fed back to z-generator 37. 

[0067] Reference is now made to Fig. 3B, which schematically illustrates the 

structure of a two-tap batch channel estimator of the present invention, here labeled 40, 

when operating on a vector y (t) whose dimension is two. Elements similar to those 

of Fig. 3 A carry similar reference numerals. 

[0068] Channel estimator 40 comprises noise variance estimator 32, a z-generator 42, 
two multipliers 38A and 38B and two summers 39A and 39B. Z-generator 42 receives 

the demodulated vector signal y(t), the probability vector p(t) from probability 

generator 34, the noise variance C from noise variance estimator 32 and the channel 
estimates h \ of the previous iteration. 
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[0069] Z-generator 42 then generates the value z(t5h') from Equation 15. 

Multipliers 38A and 38B multiply the outputs yg(t)and y^t), respectively, with 

the value z(t;h') and summers 39A and 39B sum the outputs of their respective 
multipliers 38A and 38B over a period of length T. Each summer 39 produces its 
updated channel estimate ll , which is also fed back to z-generator 42. 

[0070] It will be appreciated that the present invention is also operative for more than 

two taps. The structure of the channel estimator is similar to that of Fig. 3B. 

[0071] The channel estimator of the present invention can also be implemented in a 

sequential (or adaptive) manner. For these implementations, the channel estimator 

updates die estimate one sample at a time. This amounts to finding a sequential solution 

to Equation 16 and it typically takes some time to converge. Adaptive solutions 

inherently assume a slowly time-varying channel so that channel variations can be 

tracked. Initialization for these implementations is as described hereinabove. 

[0072] Reference is now made to Fig. 4, which illustrates a general adaptive channel 

estimator 100. Estimator 100 is similar to estimator 30 of Fig. 2 and comprises an 

adaptive channel tap determiner 102 in place of channel tap determiner 36. 

[0073] Adaptive channel tap determiner 102 sequentially solves Equation 27 

(hereinbelow) where, for each new time instance t, a new solution is obtained. 



h, = iy(m)z(m;h,)*A, w 

m=l 

Equation 27 
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X is an exponential forgetting factor. With the special choice of X = 1 , 
Equation 27 reduces to Equation 16 up to a scalar gain. 

[0074] There are a variety of sequential algorithms, two of which are presented 
hereinbelow. 

[0075] A first approach uses anchors and linear interpolation between adjacent 
anchors similar to the one described in the following documents that are incorporated 
herein by reference: 

[0076] H. Andoh, M. Sawahashi, and F. Adachi, "Channel Estimation Using Time 
Multiplexed Pilot Symbols for Coherent Rake Combining for DS-CDMA Mobile 
Radio," Proceedings of the IEEE Vehicular Technology Conference, 1997, pp. 954 - 
958, 

[0077] F. Adachi and M. Sawahashi, "Wideband Wireless Access Based on DS- 
CDMA," IEEE Transa ctions on Comm unications, pp. 1305-1316, July 1998. 
[0078] R Adachi, M. Sawahashi, and H. Suda, "Wideband DS-CDMA for Next 
Generation Mobile Communications Systems," IEEE Communicat ions Mazarine, 
September 1998. 

[0079] In the prior art anchor approach, illustrated in Fig. 5 to which reference is now 
briefly made, the channel estimator averages a few pilot symbols from the beginning of 
each of a series of successive time-slots to generate a sequence of "anchors" 50 and 52. 
Then, the channel estimator linearly interpolates between the anchor at the beg innin g of 
one slot to that at the beginning of the next slot to generate the channel estimates 54. 
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[0080] In accordance with a preferred embodiment of the present invention, the channel 

estimator computes a pilot anchor h p using the N p pilot symbols of one time slot and 
a data anchor h s using N g data symbols, as follows: 

h p =-^--Ey(k).s(k)* 

2N p k=i 

Equation 28 

h s =^-Zy(k)-z(k;h p )* 

2N S k=i 

Equation 29 

where N. is a user-defined parameter determining the number of data 
symbols to be used One may set N s = 2N p so that N p of the data symbols are 

taken from before the pilot symbols and N p of the data symbols are taken from after 
the pilot symbols. The channel estimator then averages the two values: 

h n +h, 
h = p s 

anchor ^ 
Equation 30 

[0081] The channel estimator linearly interpolates between adjacent anchors 

h anchor ( n ~~ -0 311(1 ^ anchor ( n ) to obtain ^ channel estimates for the n -th slot. 
In accordance with a preferred embodiment of the present invention, the channel 
estimator separately interpolates the amplitudes and phases of the adjacent anchors. 
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[0082] In order to improve upon the statistical variability of the anchors, the older 

anchor can be averaged, as follows: 

M ^ 

£h mchor (n-l-k)-p k 

H anchor ( n ~1) - — ~ m 

k=-M 

Equation 31 

and, as before, phase and amplitude linear interpolations between 
^■anchor C 11 ~~ *0 311(1 ^anchor ( n ) generate the II -th slot channel estimator. Both 
M and the u k 's are user defined parameters. For example, the values M = l. 

= fi l =0-8 and J3 Q = 1 were used in H. Andoh, M. Sawahashi, and F. Adachi, 

"Channel Estimation Using Time Multiplexed Pilot Symbols for Coherent Rake 
Combining for DS-CDMA Mobile Radio, 3 ' Proceedings of the IEEE Vehicular 
Technology Conference, 1997, pp. 954 - 958. 

[0083] Reference is now made to Fig. 6, which illustrates another adaptive channel 
estimator, labeled 60. Channel estimator 60 is similar to channel estimator 30 of Fig. 
3 A except that summer 39 is replaced by an infinite impulse response (HR) filter 62. 

[0084] In this embodiment, the product y (t) • z(t \ h) is filtered with some pre- 
specified IIR filter. Here, we have the channel estimator given by: 

h(t)--h(t-1)a 1 -...-h(t-p)a p +y(t).z(t;h^ 
+y(t-1)z<t-1;h(t-2))^ 
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Equation 32 

where z(t^h(t — 1)) is given by Equation 17 for pilot symbols, and by 
Equation 19 for data symbols (or by a similar expression for TPC bits). The filter 
coefficients 3,j and bj are user-defined parameters as are the filter orders pand q. 
Note that taking p = q = 1 reduces to the well known exponential forgetting factor 
sequential solution of Equation 32 where the iterations are performed sequentially in 
time. In general, the specific choice of filter parameters, should represent the best 
compromise between noise rejection and filter delay. In order to maintain a low 
complexity implementation, the data should be decimated prior to applying Equation 32. 
For example, for a constant input rate of 16ksps, one of 4 samples is used in Equation 32 
when the symbol rate is 64ksps. However, the pilot symbols should not be decimated. 
Thus, for a slot containing 4 pilot symbols and 36 data symbols, the estimator uses 10 
samples per slot, 4 of which correspond to pilot symbols and the other 6 correspond to 
data symbols. The data can also be averaged prior to the filtration operation, rather than 
decimated. A further alternative is to combine the two approaches and use averaged 

symbols only at the feedback part of Equation 32 (i.e. the terms multiplied by the 2L { ). 

This produces a very slowly varying, adaptive channel estimator but does not increase 
the estimator's delay. Such slowly varying estimators can be utilized to track slowly 
fading channels. 

[0085] The methods and apparatus disclosed herein have been described without 
reference to specific hardware or software. Rather, the methods and apparatus have 
been described in a manner sufficient to enable persons of ordinary skill in the art to 
readily adapt commercially available hardware and software as may be needed to reduce 
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any of the embodiments of the present invention to practice without undue 

experimentation and using conventional techniques. 

[0086] It will be appreciated by persons skilled in the art that the present invention is 
not limited by what has been particularly shown and described herein above. Rather the 
scope of the invention is defined by the claims that follow after the appendix: 
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APPENDIX A 

[0087] There are several simple approximations that can be used in Equation 19 in 
order to simplify the implementation. First we may replace the operands of the sinh 
operator as follows: 

\jteal{<m>h*>Q-j)} 
a 



4-MaxJRea/{^ 



-\ReaI{<y(t)X>(} + M 
a 

»-^A£axJpcalfa(t>K ■&+/)]• Signal fa(t) -h t * -(l+j)}] 

Equation A-l 

[0088] Note that this approximation simply uses the finger with the highest temporal 
SNR to estimate the transmitted symbol. This is a well-known selection criterion used in 
numerous diversity combining scenarios and contrasts with the original criterion of 
Equation 19, which is a maximal ratio diversity combination criterion. If we further 
approximate and consider only the maximum of the two terms in Equation A-l 9 then 
Equation 1 9 reduces to 

z(t;h ML )«s(t).tanh(x(t)) 
Equation A-2 

where 

Equation A-3 
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[0089] A further approximation can be implemented when considering the case of 

unequal noise variances, i.e. the case where each component of n(t) possesses a 

different variance. The log-likelihood function of Equation 10 should, in this case, be 

replaced by 

-T'ly.W-h.-d+jJr/a 1 , 



P(y(t);h)=c 



Pi(t)-e 



Dim (7) , 

£ |y k (t)-vo-j)| 
-Y'ly.W-M-i-.or/a 1 . 



+ p 2 (t)-e 
+ p 3 (t).e 

+ P 4 (t)-e 



Equation A-4 
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